Critical sets of the total variance of state detect all SLOCC entanglement classes 
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We present a general algorithm for finding all classes of pure multiparticle states equivalent 
under Stochastic Local Operations and Classsical Communication (SLOCC). We parametrize all 
SLOCC classes by the critical sets of the total variance function. Our method works for arbitrary 
systems of distinguishable and indistinguishable particles. We also show how to calculate the Morse 
indices of critical points which have the interpretation of the number of independent non-local 
perturbations increasing the variance and hence entanglement of a state. We illustrate our method 
by two examples. 
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The problem of classifying pure states up to Stochas- 
tic Local Operations assisted by Classical Communica- 
tion (SLOCC) has been intensely studied during the last 
decade by many authors jH-Q. Although the SLOCC 
classes are known for some particular systems, for ex- 
ample in the cases of three or four qubits, the general 
method allowing similar derivations for arbitrary system 
of many particles which treats in the unified way distin- 
guishable and indistinguishable particles has been miss- 
ing. In this letter we provide such an algorithm. It is 
based on studying the structure of critical points of the 
total variance of a state with respect to local unitary (LU) 
operations. Since the total variance of the state can be 
interpreted as an entanglement measure Q our approach 
has a clear physical meaning. Although our results are 
based on some relatively advanced mathematical tools, 
the final algorithm we propose reduces the problem to 
simple computations involving diagonalized reduced den- 
sity matrices. 

The paper is organized as follows. First we give a defi- 
nition of the total variance of a pure state, Var( [ip] ) . Next 
we discuss a connection between the concept of the mo- 
mentum map and the reduced one-particle density ma- 
trices. Following this we show that the total variance 
of a state is, up to some unimportant additive constant, 
the norm squared of the momentum map. This observa- 
tion allows us to identify the critical points of Var ([ip]), 
i.e. points [i[>] in F(H) for which dVax([ip}) = 0, with the 
critical points of the momentum map. The later turn 
out to be well understood in Geometric Invariant The- 
ory (GIT) d, 0j- Using its results we show that for each 
SLOCC class of pure states the restriction of Var ([?/>]) 
to it attains maximum on exactly one orbit of the local 
unitary action. This orbit contains the most entangled 
representatives of a given SLOCC class. Finally, we de- 
scribe how to calculate the Morse index of Var([^)]) at 
critical points and explain that it has a meaningful inter- 
pretation as the number of independent non-local pertur- 
bations which increase the total variance, Var( [ip]) thus, 
consequently, also its entanglement. The whole idea is 
illustrated by familiar examples. 



Throughout the paper we consider the system of 
L identical particles which can be either distinguish- 
able or indistinguishable. For the former the Hilbert 
space is the tensor product of L copies of C^, Hd = 
C N !g) . . . £g> C N . Local unitary operations are repre- 
sented by the direct product of L copies of SU(N), 
K = SU(N) xL , and SLOCC operations by of L copies 
of SL(N), G = SL(N) xL , (see (J^). The action of 
an element (Ai, . . . , At,) of either K or G on the vector 
&T-Ld is given by 



I—?- A\ ® . . .®A L \ip) 



(1) 



For indistinguishable particles, i.e. for systems of 
L bosons or L fermions the relevant Hilbert spaces 
are fully symmetric or fully antisymmetric parts of 
the full tensor product, i.e. Hb = Sym L (C N ) and 
Hf = f\ L (C w ). The local groups corresponding to 
LU and SLOCC operations are K = SU(N) and 
G = SL(N) and the action of an element A E K (or 
A e G) on \tp) G H b /f is given by ((TJ with A k = A. 
Because the global phase factor and the normalization 
are physically irrelevant it is useful to identify pure 
states with elements of the complex projective space 
F(H). To this end we identify vectors from H that 
differ by a nonzero complex scalar factor and denote by 
[ip] £ F(H) the state corresponding to vectors G H 
under this identification. Because groups K and G act 
linearly on T~L they also act naturally on P('H). Two 
pure states [ip] , [0] G P (H) are SLOCC equivalent if 
and only if \ip) and \<p) can be connected by the action 
of G on the level of H, i.e. g.\ip) = e iQ |0), a G [0, 2tt), 
for some g G G. Actions of K and G on H induce 
actions of the corresponding Lie algebras t and g. 
The action of a = (ai, . . . , cx-l) in either su(N)® L or 
sl(N)® L on the vector G H is given by a.\ip) = 
(ai ® I N ® ■ • ■ ® In- + ■ ■ ■ + In ® • ■ • <£> In ® ot L ) 
For indistinguishable particles a k = a G su(A^) (or 
a G sl{N)). 

After fixing notation we can define the total variance 
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of the state [ip] € V(H), first introduced in [j| 



Hilbcrt-Schmidt norm. It is easy to see that: 



dim K 



Vax(M) 



dim A" 



(2) 



where {X^f"™ is an orthonormal (with respect to the 
Hilbert-Schmidt scalar product (A, B) = Tr(A* B)) basis 
of t represented as matrices acting in % [IH . The expres- 
sion (|2|) is simply the sum of variances of 'local' observ- 
ables Xi calculated in the state [ip]. It can be checked 
that Var([ , 0]) is K - invariant and attains its minimum 
precisely on the set of separable states. In fact, it can be 
used as a measure of entanglement 

In order to proceed we briefly describe the concept 
of the momentum map. This is a map which encodes 
information about all first integrals of a given classical 
Hamiltonian system with symmetries. For example, for 
a Hamiltonian system with the rotational symmetry, i.e. 
when the symmetry group is K = 50(3), the first in- 
tegrals are three components of the angular momentum 
corresponding to the invariance of the system with re- 
spect to infinitesimal rotations along axes x, y, z. The 
infinitesimal rotations generate the Lie algebra so (3) of 
50(3). There are many possible choices of basis for so(3), 
each corresponds to different choice of rotation axes and 
each gives three first integrals. The mathematical ob- 
ject which encodes information about all first integrals 
for all possible choices of generators of so(3) is a map 
(i : M — > so(3)*, i.e. a map from the phase space M to 
the space of linear functionals on the Lie algebra 50 (3) 
of the group 50(3). For every infinitesimal symmetry 
£ € so (3) one can find the corresponding first integral by 
evaluating fi^(x) = /j,(x)(£). This idea can be generalized 
to arbitrary K and the corresponding map fx : M 6* 
is called momentum map 0. 

It is well known that P (H) posses a natural phase space 
structure [19j . Moreover, for considered composite sys- 
tems we can treat LU operations as symmetries of the 
system as they do not change entanglement. Remark- 
ably, the momentum map for both distinguishable and 
indistinguishable particles has a clear quantum mechan- 
ical meaning. It is a map which assigns to each state its 
reduced one-particle density matrices. For distinguish- 
able particles we have 

M) = - jjlN, ■ • • . PirW) - ^) . (3) 

where Pi([ip]) is the normalized i-th reduced one-particle 
density matrix of the vector state | ip). For indistinguish- 
able particles we have: 



MM) = piM) - jjIn 



(4) 



IIMM)ll a = 



(5) 



i=i 



We also notice that C2 = J2i=i K ^ s the representa- 
tion of the second order Casimir operator [i| , a general- 
ization to other Lie algebras of the squared total angular 
momentum relevant for the above mentioned example of 
rotationally invariant system. The operator C2 commutes 
with each Xi. Since the group K acts on H irrcducibly, 
it is proportional to the identity operator with the pro- 

(ip\e \ib) 

portionality constant c = (\p^p) ■ Making use of this 
formula and ([5]) we get the final expression for the total 
variance of a state: 



Var(M) = c-|| M (M)|| 



(6) 



Since fJ-([ip]) is an element of t we can calculate its 



Notice that by @ the critical points Var( [?/>]) are exactly 
the critical points of 1 1 /i ( [ip] ) 1 1 2 . The later where intensely 
studied in GIT context in the 80'. Specifically, the fol- 
lowing theorem is the reformulation of the general result 

Theorem 1 Assume that [ip] is a critical point of Var : 
P('H) — > M. Then the restriction of the total variance 
function Var to the SLOCC class of [ip], i.e. G-orbit 
through [ip], attains its maximum value on a unique K- 
orbit, which is the orbit K.[ip] through the state [ip]. 

In general it may happen that a particular G-orbit does 
not contain any critical point of Var : ¥(%) — s- M. Never- 
theless, as we show in [101 ] . such an orbit always contains 
in its closure the if-orbit which is critical for Var ([-(/>]). 
The mathematical details of this construction are rather 
subtle and wc will not discuss them here. We would like 
to emphasize, however, that knowing critical if -orbits 
allows the reconstruction of these G-orbits which do not 
contain any critical point of Var (see [l(|). Under the 
above assumptions Theorem [T] can be interpreted as the 
statement saying that the critical sets of the total vari- 
ance function of a many-particle system parametrize all 
SLOCC classes of states. Moreover the distinguished 
local unitary orbit on which Var attains a maximum 
contains maximally entangled representatives of a given 
SLOCC class. Although we already know that the to- 
tal variance cannot increase when we perturb the state 
from this orbit with SLOCC operations, it is still not 
clear what happens when the perturbation is non-local. 
This information is stored in the Morse index of a crit- 
ical point i.e. the number of negative eigenvalues of the 
Hessian of \\fi ([ip]) || 2 calculated at [ip]. We will investi- 
gate this idea on two examples. First, however we need 
some effective way for finding critical points of ||^([-0])|| 2 
and calculating the Morse index. In [lj| we prove the 
following theorem 
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Theorem 2 1. The state [ip] is a critical point of 
Var : P('H) — > R if and only if the corresponding 
vector [ip) is an eigenvector of p,([ip]), i.e. for some 
X € R we have 



p([ip])]ip) = X]ip). 



I?. The Morse index of \\p 



(7) 

at a critical point [ip] 



is the the Morse index of the function f : % — > 



(8) 



calculated at \(p) = \ip) and restricted to the orthog- 
onal complement (with respect to the usual inner 
product onH) of the subspace Q.\ip) 

In Equations © and §5$ p([ip]) is understood as an op- 
erator acting on H. Expression Q.\ip) denotes the image 
of the vector \ifj) under the action of the whole g rep- 
resented on H, where g — sl(N)® L for distinguishable 
particles and q ~ sl(N) for indistinguishable particles. 

Using ([7]) we can divide the problem of finding critical 
points of Var( [ip]) into two conceptually different parts. 
First, we notice that condition ([Jj is automatically satis- 
fied if /u([V>]) = 0. This corresponds to states with maxi- 
mally mixed reduced one-particle density matrices. Since 
the set ^t -1 (0) is if-invariant, each if -orbit contained in 
^i" 1 (0) represents a different SLOCC class of states. By 
splitting p~ x (0) into disjoint orbits of K and picking ex- 
actly one state from each orbit we get a parametrization 
of all SLOCC classes of this type. These classes con- 
tain almost all states [l(J. For effective calculations it 
is desirable to have some canonical form of the state. 
In iV-particle case such a form was given in The 
remaining critical points of Var(V>) stem from ([7]) when 
A t ([V']) 7^ 0- The following two remarks are crucial in this 
case (see [HI for a detailed discussion) 

1. For any \tp) there exists LU operator U = U\ ® 
. . . <X> Ul such that p([ip']), where \ip') = U\ip), is given 
by (O with diagonal reduced density matrices Pi([ip']) 
whose eigenvalues are arranged in non-increasing order; 
the matrices Ui € SU(N) diagonalize pi([ip]). Notice 
that since critical sets of Var ([ip]) are K - invariant it 
implies that we can consider equation (jTj) assuming that 
ip posses the above described property of ip'. For bosons 
and fcrmions the relevant operators U also exist and are 
given by U = U x <g> . . . <g> Ux, where Ui £ SU(N) diag- 
onalizcs px([i>]) from (j4|) and makes its spectrum non- 
increasingly ordered. 

2. Let Pi(ip) be a point in M. N whose coordinates are 
given by non-increasingly ordered spectrum of Pi([ip]) — 
jfI N . The set *(P(W)) = {(Px(iP),...,PlW) ■ M e 
P("H)} is a convex polytope in M. LN , the so-called Kirwan 
or momentum polytope. This fact is a simple corollary 
[l2^ from the general theorem known as the convexity 
property of the momentum map (l3| . Moreover, in the 



variety of settings for distinguishable and indistinguish- 
able particles the inequalities defining this polytope are 
explicitly known 3, 15^ . 

Using the above remarks the search for remaining crit- 
ical points of Var([V->]) can be reduced to the following 
procedure. For each point P = (Pi, . . . , Pl) in the poly- 
tope ^(¥(H)) \ {0} we construct corresponding operator 



ap = ap 1 



in 



in 



I 



N 



I 



N 



where otp i is a diagonal matrix with diagonal elements 
given by Pi. Next we are looking for states [ip] for which 
the condition 



a P \*l>) = \\tl>) 



(9) 



is satisfied and, at the same time, = ctp. The prob- 

lem of finding critical points is thus reduced to a study of 
eigenstates of the operator a p which is an N L x N L di- 
agonal matrix for P E ^'(P('H)). If the diagonal elements 
of ap are nondegenerate the corresponding eigenvectors 
are separable states. Hence in order to find nontrivial 
SLOCC classes we are interested in the situation when 
the spectrum of ap is degenerate. Remarkably, the di- 
mensions of the degenerate cigenspaccs arc typically rel- 
atively small and hence it is quite easy to verify for which 
states /id/0]) = ap. Finally, as we explain in [10( it is al- 
ways possible to divide critical points found in this way 
into finite number of families, where each family contains 
SLOCC nonequivalent classes of states characterized by 
fixed reduced one-particle density matrices. 

In order to demonstrate how to use the above method 
we provide calculations for two simple examples (three 
qubits and three 5-state fermions). For cases when cal- 
culation are more elaborate see [10| . 

Three qubits. For the three-qubit system the Hilbcrt 
space is H = C 2 ® C 2 <£> C 2 , the SLOCC operations are 
represented by the group G — SL(2) x3 , and the local 
unitary operations by the group K = SU(2) x3 . Using 
the cannonical form of a three-qubit state 
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\ip) =p|011) + g|101) +r|110) + s|lll) + z|000), (10) 

the reduced one-qubit density matrices can be easily cal- 
culated. It is then a matter of straightforward calcu- 
lations to see that there are exactly two states of the 
form (|10p which have maximally mixed reduced one- 
particle density matrices. One of them is the GHZ state, 
\^ghz) = (|000) + [111)), and the other one is LU 
equivalent to it. Effectively, the main SLOCC class is 
the G-orbit through GHZ. The remaining classes can 
be found by considering Equation (|9|), where ap is 8 x 8 
diagonal matrix. The Kirwan polytope is given by three 
polygonal inequalities [l3| for the smallest eigenvalues 
of the reduced one-qubit density matrices and hence is 
three-dimensional. By direct calculations we find that 
the only eigenspaccs contributing new critical points are 
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of dimension 1, 3 and 4. They give separable, biseparable 
and W SLOCC classes, respectively, i.e. 



\4>BSl) = 



\4>w 
1 



1 



(1 100) + |010) + |001>) 



V2 

\lpBS3) 



V3 

(1 100) + |111)), 



V2 



(|010) + |111)) : 



V2 



(|ooi) + |iii)), i'ssp = |ooo). 



Since the state GHZ is in /x -1 (0) the Morse index at 
it is equal to zero [loj . This means it is impossible to 
increase the total variance Var by any perturbation of 
K.\ipGHz)- The case of \ipw) is more interesting, namely 
for the orthogonal complement of Q-\il>i) we have 

( #i)) ± =Span{|lll),*|lll)}. 

It is easy to check that the Hessian of the function 
f{\i/j)) defined by Eq. is negative on (Q.\ipi)) and 
hence the index equals 2. Similar calculations show that 
for biseparable states, i.e. \ipBSi), \i>BS2) and \1pBS3} 
we have index(|^ssfc)) = 6 and for separable state 
index(IV'ssp)) = 8. 

Three five-state fermions. As a second example let us 
consider a system consisting of three fermions with five- 
dimensional single particle Hilbert spaces, TL = /\ 3 (C 5 ) , 
K = SU (5), G = SL(5, C). In C 5 we fix an orthonormal 
basis {|1), |2), |3), |4), |5)} and chose the (also orthonor- 
mal) basis of A 3 (C 5 ) 

*2, «3> := \h) A \i 2 ) A \i 3 ) , 1 < h < i 2 < 13 < 5 . (11) 

Any \tp) G % has a decomposition 



\tp)= ck u i 2t i 3 \ii, 12,13) 

l<ii<i 2 <i3<5 



(12) 



where the scalar coefficients fulfill the normalization con- 
dition. A system with TL = f\ 3 (C 5 ) can be considered 
as consisting of two holes and the structure of the spec- 
trum of the rescaled density matrix p ( [tp] ) is well known 
1 51 ] - First condition says that the eigenvalues of p are 
bounded from above by |. The second one is that the 
maximal eigenvalue X ma x equals exactly i . The last con- 
dition states that all eigenvalues except for X ma x are, at 
least, doubly degenerate. Using these properties it can 
be easily shown that there are only two noncquivalent 
critical sets of Var. They are parametrized by the states: 



|1,2,3), |</>2) 



1 



— (11,2,3) + 1 1,4, 5)) . (13) 



Using Theorem 2 we find, in the same manner as in 
the case of three qubits, that index (|^i)) = 6 and 
index (\1jj2)) = 0. The latter fact is a consequence of the 
observation that (g.j^)) -1 = 0. It follows that G.^2] is 



of the maximal dimension and that the total variance at- 
tains the global maximum on the K orbit passing through 
[ip2\- It is interesting to note that in the case considered 
/i _1 (0) is empty, that is there are no states in P('H) that 
have maximally mixed reduced density matrix. 

Summarizing, we provided an algorithm for finding and 
classifying SLOCC equivalent pure states. The algorithm 
is effective and universally applicable in various setting 
including distinguishable and indistinguishable systems 
of particles with arbitrary number of components. 
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After completion of this work, we have learned about 
independent related work by Walter, Doran, Gross, and 
Christandl [T^- 
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